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Objective: Resting-state networks (RSNs), including the default mode network (DMN), have been considered as markers of 
brain status such as consciousness, developmental change, and treatment effects. The consistency of functional 
connectivity among RSNs has not been fully explored, especially among resting-state-related independent components 
(RSICs). 

Materials and Methods: This resting-state fMRI study addressed the consistency of functional connectivity among RSICs as 
well as their spatial consistency between 'at day 1' and 'after 4 weeks' in 13 healthy volunteers. 

Results: We found that most RSICs, especially the DMN, are reproducible across time, whereas some RSICs were variable in 
either their spatial characteristics or their functional connectivity. Relatively low spatial consistency was found in the basal 
ganglia, a parietal region of left frontoparietal network, and the supplementary motor area. The functional connectivity 
between two independent components, the bilateral angular/supramarginal gyri/intraparietal lobule and bilateral middle 
temporal/occipital gyri, was decreased across time regardless of the correlation analysis method employed (Pearson's or 
partial correlation). 

Conclusion: RSICs showing variable consistency are different between spatial characteristics and functional connectivity. 
To understand the brain as a dynamic network, we recommend further investigation of both changes in the activation of 
specific regions and the modulation of functional connectivity in the brain network. 

Index terms: Resting state network; Default mode network; Independent component analysis; Consistency; Functional connectivity 
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INTRODUCTION 

The resting-state network (RSN) is a term referring to 
functional networks of brain regions that are active when 
an individual is not focused on the outside world. The RSN 
is characterized by coherent oscillations at a rate lower 
than 0.1 Hz among brain regions. Generally, the default 
mode network (DMN) comprises a set of brain regions 
including the cingulate cortex, hippocampus, and medial 
prefrontal and inferior parietal cortices and includes 
members of several RSNs (1, 2). The DMN is active during 
sleep (3, 4) and under anesthesia (5), as well as when at 
rest (6), but is not activated during passive viewing with a 
natural scene (7). The DMN can be measured using various 
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neurophysiological or functional neuroimaging tools such as 
electroencephalography (EEG) (8), magnetoencephalography 
(MEG) (9, 10), functional magnetic resonance imaging 
(fMRI), and functional near infrared (fNIR) (11). The 
activity of the DMN detected with functional neuroimaging 
is similar to that of the electric or optic signal measured 
directly from the brain cortex (12, 13). While performing a 
task, the DMN shows a reciprocal relationship with the task- 
related active network (14). This means that the activity 
of the DMN may be decreased while performing a task and 
increased during the resting state. A previous study reported 
that the more extensive the cognitive task is, the more 
active the subsequent resting state of the DMN becomes 
(15). Thus, the RSN, including the DMN, might form part of 
the backbone of brain behavior. 

It is understood that the DMN plays a role in cognitive 
processes such as self-reflection, moral reasoning, 
recollection, or imagining the future (6). The DMN may also 
be associated with some psychopathologies, (e.g., self- 
disturbance in schizophrenia) (16), and inattentiveness 
in attention deficit hyperactivity disorder (ADHD) (8, 17). 
In previous neuroimaging studies, deficits in the DMN 
have been reported in various neuropsychiatric disorders 
(18), including schizophrenia (10, 14), autism (19), ADHD 
(20-23), Alzheimer's-type dementia, and mild cognitive 
impairment (24-26). From a clinical perspective, a question 
of interest is whether a deficit in the DMN can be restored 
or modulated. Prior studies have reported that deficits in 
the DMN can be modulated with antipsychotic medication in 
schizophrenia (27) and with psychostimulants in ADHD (28, 
29). These results suggest that the characteristics of the 
DMN may provide a measure of treatment effectiveness as 
well as confirmation of diagnosis. To measure the treatment 
effect, however, the characteristics of the DMN need to 
be consistent across time. Previous studies using the fMRI 
data of healthy individuals have reported the test-retest 
reproducibility of the spatial characteristics of the RSN 
including the DMN (30, 31). However, the reproducibility of 
the functional connectivity of the RSN, including the DMN, 
has not yet been established using fMRI data, even though 
the functional connectivity among resting-state-related 
independent components (RSICs) has been used to explore 
clinical significance (32). 

Usually, the DMN is constructed using a correlation map 
between the time course of a seeding point or region 
of interest (ROI) and those of all other voxels within 
the brain (33, 34). This seed-based correlation analysis 
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using anatomical ROIs could suffer from the confounding 
factor of individual anatomical variability (35), and could 
also show only relationships of one-to-multiple ROIs and 
not relationships among all ROIs. An alternative ways 
of constructing the DMN is to extract raw time-course 
data from the ROI selected on the spatial map of the 
independent component (IC) (32). However, the raw time 
courses might also be affected by various forms of noise, 
even though the seeding point was chosen based on the IC 
map. In the present study, the time courses of RSICs, rather 
than raw time courses, were used for measuring functional 
connectivity. Our approach might be less affected by either 
individual anatomical variability, as with a seed-based 
approach, or by noise, as with an ICA approach using raw 
time courses. Also, all connections among RSICs could be 
explored with our approach using the time courses of RSICs, 
unlike the seed-based correlation analysis. In this resting- 
state fMRI (rs-fMRI) study, we addressed the consistency 
of functional connectivity among the time courses of RSICs 
and also confirmed the spatial characteristics of RSICs in 
healthy individuals. 

MATERIALS AND METHODS 

Subjects 

Enrollment for subjects was advertised on the message 
board at Eulji University Hospital, Daejeon, Korea. Thirteen 
healthy young male subjects (Age = 22.4 + 2.5 years) 
underwent fMRI. Subjects were interviewed by a psychiatrist 
using DSM-IV criteria to exclude any possible psychiatric 
disorders. All subjects were right-handed and had no history 
of psychiatric or neurological disorders or head trauma. This 
study was approved by the institutional review board of 
the Eulji University Hospital. All subjects provided written 
informed consent before participation in the present study. 

Functional MRI Experiments 

Subjects underwent fMRI scan of 5 minutes and 30 
seconds both on day one and after 4 weeks at rest with 
their eyes closed. Using 3.0 Tesla MRI (IS0L, Gwangju, 
Korea), fMRI was acquired with a BOLD-sensitive echo- 
planar gradient-echo (EPI) sequence with the following 
imaging parameters: repetition time (TR), 2000 ms; echo 
time (TE), 35 ms; flip angle (FA), 70°; spatial resolution, 
3.4375 x 3.4375 x 5 mm 3 , imaging matrix, 64 x 64; field- 
of-view (F0V), 220 x 220 mm 2 , number of slices, 25, 165 
volumes. Subjects were instructed to keep their eyes closed 
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and not to think of anything particular during fMR scans. MELODIC (Multivariate Exploratory Linear Optimized 

Decomposition into Independent Components) module of 
Data Preprocessing the FSL software (http://www.fmrib.ox.ac.uk/fsl/feat5/ 

The fMRI data were processed and analyzed using the index). The first five scans were discarded to account for Tl 
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Fig. 1. 22 RSIC maps and their highest spatial correlation coefficients (C.C.) among possible 63 IC pairs between two time points 

(at day 1 and for after 4 week). Posterior default mode network (IC 1), visual system (IC 2, IC 3, IC 10, IC 15, IC 17), sensorimotor system (IC 4, 
IC 5, IC 6), subcortex (IC 7, IC 22), frontoparietal network (IC 8, IC 21), superior parietal network (IC 9), anterior default mode network (IC 11, 
IC 20), temporal network (IC 12), executive control system (IC 13, IC 14, IC 16), cerebellum (IC 18), motor system (IC 19), parietal network (IC 
21). 
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saturation effects. The remaining 160 images were spatially 
realigned using a rigid-body transformation and were 
subjected to slice-timing correction. Next, a brain mask 
from the first fMR data volume was created to eliminate 
signals outside of the brain in each subject. Spatial 
smoothing using a full-width Gaussian kernel at half- 
maximum, 6 mm was performed to reduce noise without 
reducing valid activation. The smoothed images were 
filtered with a 128-second high pass filter to effectively 
remove signal drift. The serial correlations were removed to 
validate the statistics and be maximally efficient. Finally, 
the functional images were registered to the Montreal 
Neurological Institute (MNI) Tl template using FLIRT 
(FMRIB's Linear Image Registration Tool) for FSL. 

Spatial Consistency of RSIC Maps Across Time 

To test the spatial consistency, multisession temporal 
concatenation (36) was run on all 13 participants as a 
group at each time point where a standard (space x time) 
ICA decomposition was conducted. Both the principle- 
component analysis (PCA) and the subsequent ICA 
algorithm produced 54 ICs at week 0 and 63 ICs at week 4. 
Each IC map was divided by the standard deviation of the 
residual noise and was plotted onto a histogram of intensity 
values (37). MELODIC was used to carry out inference 
on the estimated maps using a Gaussian/gamma mixture 
model and an alternative hypothesis testing approach at 
a threshold level of 0.5. A threshold level of 0.5 in the 
case of alternative hypothesis testing means that a voxel 
survives thresholding as soon as the probability of it being 
in the active class exceeds the probability of it being in the 
background-noise class (37). The spatial cross-correlation 
analyses among 3402 IC pairs (54 x 63), between two time 
points, were performed with thresholded IC maps using FSL. 
This means that 63 correlation coefficients were produced 
for each of the 54 ICs on day 1. The IC pair showing the 
highest correlation coefficient among the 63 possible ICs 
was considered as the same IC between two time points 
and their correlation coefficients used for a value of spatial 
consistency. Among the 63 ICs, 41 were not RSICs were 
signals from white matter, cerebrospinal fluid, vessels, and 
motion artifacts. The remaining 22 pairs were selected as 
the RSIC and were ordered with their correlation coefficients 
(Fig. 1). The spatial correlation coefficients were compared 
between RSICs and non-RSICs using a p value < 0.05 as a 
significance threshold. While this approach provided the 
degree of spatial consistency of each IC, it did not provide 



Jeong et a I. 

its statistical difference between two time points. 

Another method using a post-hoc regression analysis 
was employed to test the question about whether spatial 
maps of RSICs are consistent across time. Unlike the first 
ICA method performed for each time point individually, this 
multisession temporal concatenation (37) was run on both 
time points of all the 13 participant pairs as a group. Both 
the PCA and subsequent ICA algorithm yielded 48 spatial 
IC maps. Each IC map was thresholded by fitting a mixture 
model using the same method outlined above (37). Among 
48 ICs, 26 that were not RSICs were discarded as signals 
from white matter, cerebrospinal fluid, vessels, and motion 
artifacts. The remaining 22 RSICs were selected for a 
pairwise group comparison (Fig. 2). For the pairwise group 
comparison, to identify differences in spatial characteristics 
of RSICs according to time point, we specified a temporal- 
design matrix, a subject-design matrix, and corresponding 
contrast matrices. Finally, a post-hoc regression analysis 
was performed on estimated time courses and subject 
modes at a statistical threshold level of p < 0.05. 

The Consistency of Functional Connectivity Among Time 
Courses of RSICs Across Time 

To extract the RSIC time courses, we used the 22 ICs 
selected in the prior step for spatial consistency using 
the Gaussian/gamma model. To identify subject-specific 
temporal dynamics within each individual's fMRI dataset, 
the full set of group-ICA spatial maps was used in a linear 
model fit against the separate fMRI data sets (38). This 
process produced time courses for each component and 
subject. Time courses of the 22 RSICs of each subject 
were selected for further functional connectivity analysis. 
Here, both cross-correlation (CC) and partial-correlation 
(PC) coefficients were used as measures of functional 
connectivity. These correlation analyses produced CC and 
PC matrices of 22 by 22 in each subject. The CC referred to 
the covariance between time courses of members in each 
IC pair. The PC referred to the correlation of time courses 
between two ICs after accounting for the relationship of 
each time course to the other 20 reference time courses. 

To test the consistency of functional connectivity across 
time, two approaches were applied. First, using correlation 
coefficients for all IC pairs, the intraclass correlation 
coefficient (ICC) was calculated between the two time 
points. The correlation coefficient of each IC pair was 
averaged across all subjects. Then, 462 pairs excluding 22 
self-pairs in the group-averaged coefficient matrix (22 by 



268 



Korean J Radiol 13(3), May/Jun 2012 



kjronline.org 



Functional and SpatiaL Consistency of Resting State- Related Independent Components in fMRI 



Korean Journal of Radiology 

KIR 



IC1 IC2 IC3 IC4 IC5 IC6 IC7 IC8 IC9 IC 10 

IC 1 1 IC 12 IC 13 IC 14 IC 15 IC 16 IC 1 7 IC 18 IC 19 IC20 IC 21 IC22 



□ ii ■ 

■ p B | — 

■ ■■■■ □ ■■ 

■ ■■■ 

II ■ 

■ ■■■ ■■■■■■ 
!■ 

mm , , 

■ ■■ ■ □ □ 

■ret \r • 

mam m 

□ 

■ I"- - ■ .\ 

rmr 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 



22 
21 
20 
19 
18 
17 
16 
15 
14 
13 
12 
11 
10 
9 
8 
7 
6 



































□ 






■ 




■ 














































■ 






■ 








































■ 














■ 
























■ 










■ 




































■ 




■ 






1 














■ 














































































■ 


















J 






■ 






























































































































































■ 














■ 






■ 












■ 
















































































■ 


















■ 


■ 






















































■ 
























■ 


















■ 








■ 












■ 


■ 




































■ 


■ 






























1 






■ 








■ 






































■ 








■ 


















































1 2 


3 


4 


5 


6 


7 8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 2 


1 22 



B C 

Fig. 2. Resting state related independent components (RSICs) and their functional connectivity maps. 

RSICs (A). Functional connectivity maps for cross correlation (B) or partiaL correlation (C) among RSICs. In each diagonal matrix, left upper part 
of matrix represents functional connectivity map at day 1, while right lower part represents functional connectivity map after 4 weeks. Black-line 
box represents independent component pair, which showed significantly higher mean correlation coefficient in corresponding time point (e.g., 
black-line box in left upper part means that correlation coefficient is higher at day 1 than after 4 weeks). 



22) were used to calculate the ICC. The Psych statistical 
package "R" (http://www.personality-project.Org/r/html/ 
ICC.html), was used to calculate the ICC between time 
points identified as "at day 1" and "after 4 weeks." The 
second approach was the comparison of the correlation 
coefficient of each IC pair between two time points. Fisher's 
r-to-z transformation was performed on the correlation 
coefficients. The z-transformed mean correlation coefficient 
of each IC pair was compared between "at day 1" and "after 
4 weeks" using a paired t test at a statistical threshold level 
of p < 0.05. Corrections for multiple comparisons and false- 
discovery rate (FDR) were applied to the p values resulting 



from the paired t test. 
RESULTS 

Spatial Consistency 

As a result of spatial cross correlation using FSL, 22 IC 
pairs showing the highest correlation coefficient between 
two time points were identified as RSICs. The spatial 
correlation coefficients among 22 RSICs ranged from 0.20 to 
0.74 (0.53 + 0.13) across RSICs (Fig. 1). The most spatially 
consistent RSIC was in the posterior cingulate cortex 
(PCC)/precuneus region known as a posterior DMN (39). 
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Relatively high spatial consistency was also found in the bilateral precentral/postcentral gyri (0.68). RSICs showing 
bilateral inferior occipital gyri (0.69), bilateral calcarine relatively low spatial consistency were found in the basal 
gyri (posterior medial visual component, 0.69), and ganglia (0.20), a parietal region of the left frontoparietal 



Table 1. Spatial Consistency of Resting State Related Independent Components (RSICs) Using Results of Post-Hoc 
Regression Analysis 
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Note. — N = network, S = system, DMN = default mode network, SMA = supplementary motor area, Inf = inferior, Sup = 
superior, C = cortex, G = gyrus, L = left, R = right, B = bilateral 
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network (0.36), and the supplementary motor area (SMA) 
(0.38). The remaining 26 non-RSICs including physiological 
noise, non-gray matter signals, and motion artifacts showed 
significantly lower correlation coefficients (0.38 + 0.17, t 
test, t = 3.58, df = 51.65, p = 0.001) than did the RSICs. 
The spatial correlation coefficients and the map of each 
RSIC are presented in the Figure 1. 

In the post-hoc regression analysis with a mixed 
Gaussian/gamma model, the spatial characteristics of each 
RSIC pair were not significantly different between day 1 and 
after 4 weeks (z = 0.10-1.21, p = 0.10-0.90) (Table 1). 

The Consistency of Functional Connectivity: ICC Across 
Time Interval 

Intraclass correlation analysis showed that the correlation 
coefficients among all pairs of IC time courses were 
significantly consistent across the time interval. The ICC of 
CC was 0.78 (F = 7.93, p < 0.001, CI = 0.74-0.81), while the 
ICC of PC was 0.75 (F = 7.05, p < 0.001, CI= 0.71-0.79). 

The Consistency of Functional Connectivity: Comparison 
of Mean Correlation Coefficients Across a Time Interval 

(Fig. 2) 

CCs between several IC pairs were significantly different 
across a time interval. The CC was significantly decreased 
between IC3 (upper part of the precuneus) and IC4 (superior 
medial frontal gyrus: anterior DMN; t = 2.08, p = 0.05), 
IC9 (bilateral middle temporal gyri; t = 2.73, p = 0.01) or 
IC19 (bilateral middle temporal/lateral middle occipital 
gyri: V5; t = 2.13, p = 0.04). A decreased CC was also 
found for the IC12 (right angular and supramarginal gyri/ 
intraparietal lobule) - IC19 pair (t = -2.20, p = 0.04). CC 
also significantly increased between IC17 (bilateral inferior 
occipital gyri) and IC6 (bilateral frontal poles/middle frontal 
gyri; t = -2.15, p = 0.04), IC7 (left inferior frontal gyrus; t = 
-2.23, p = 0.04), IC8 (bilateral supramarginal gyri; t =-2.48, 
p = 0.02), and IC12 (t = -2.29, p = 0.03). Increased CCs 
were also found for pairs IC9-IC21 (left precentral gyrus) (t 
= 2.49, p = 0.01) and IC15 (bilateral calarine gyri: anterior 
medial visual component) - IC20 (bilateral precentral/ 
postcentral gyri) (t = 2.15, p = 0.04). In terms of PCs, 
the functional connectivity was significantly decreased in 
the IC12 - IC19 pair (t = -2.13, p = 0.04) and increased in 
the IC18 (bilateral cuneus/calcarine gyri: posterior medial 
visual component) - IC22 (right precentral/postcentral 
gyrus) pair (t = 2.63, p = 0.01). The functional connectivity, 
both for CC and PC, of the precuneus/posterior cingulate 
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cortex component (IC 1, posterior DMN) with other RSICs, 
did not significantly change over 4 weeks. The differences 
in CC or PC disappeared after the correction for multiple 
comparisons using FDR. 

DISCUSSION 

This rs-fMRI study evaluated the consistency of both 
functional connectivity and spatial characteristics 
between 'at day 1' and 'after 4 weeks' among RSICs. First, 
we confirmed the spatial consistency of RSICs (30, 31) 
using two different methods: correlation analysis and 
the Gaussian/gamma model. Also, we found that spatial 
consistency was variable across RSICs. Second, our high 
intraclass correlation coefficient (ICC) result suggested 
that the overall functional connectivity among RSICs is 
consistent across time. However, in the comparisons of 
the correlation coefficients of each IC pair across time, the 
functional connectivity of several RSIC pairs was variable 
across time. Finally, RSIC pairs showing differences were 
more frequent in the cross-correlation (CC) analysis using 
Pearson's coefficient than in the PC analysis. Thus, our 
results suggest that both the spatial map and functional 
connectivity are consistent across time, but that the degree 
of their consistencies is variable across RSICs or by the 
correlation analysis method. 

Our results regarding the spatial consistency of RSICs 
are consistent with previous results. We should also note 
that the spatial consistency was variable across RSICs (40), 
although they were more spatially consistent than were 
non-RSICs. Relatively high consistency was found in the 
posterior DMN, bilateral inferior occipital gyri, bilateral 
calcarine gyri (posterior medial visual component), and 
bilateral precentral/postcentral gyri, whereas relatively low 
consistency was found in the basal ganglia, parietal region 
of the left frontoparietal network, and supplementary 
motor area (SMA) (Fig. 1). A previous study using ICA 
also reported both variability of spatial consistency and 
relatively high spatial consistency of RSICs compared with 
non-RSICs (31). In that prior study, RSICs showing relatively 
high consistency were the bilateral frontoparietal network, 
bilateral occipital poles (bilateral inferior occipital gyri), 
posterior DMN, bilateral inferior frontal gyri, and bilateral 
supramarginal gyri, whereas the cerebellum, SMA, and basal 
ganglia had low consistency. The spatial consistency of 
most RSICs in the present study was similar to this prior 
study, in which the sessions were 45 min or 11 months 
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apart. Both the present and the prior studies suggest that 
DMNs are spatially reproducible and that the low spatial 
consistency of RSICs in both SMA and basal ganglia should 
be further investigated in a longitudinal study. The low 
spatial consistency of the frontoparietal network, especially 
in the left parietal region, was not consistent with the prior 
study. However, there has been an rs-fMRI study supporting 
our results (41), and they reported a systematic impairment 
of associative frontoparieto-cingulate areas in altered states 
of consciousness. The discrepancy between these studies 
regarding spatial consistency in the frontoparietal network, 
especially in the left parietal region, should be investigated 
using a high-dimensional ICA analysis, which could separate 
the frontoparietal network into frontal and parietal regions. 

We also note the variability of functional connectivity 
in several RSICs, as well as the overall consistency of 
RSICs across time. In our study, the RSIC pairs showing 
differences in functional connectivity were mainly the 
connections within visual RSICs or between visual and 
other RSICs (Fig. 2). A MEG study has also shown that 
graph metrics of functional connectivity were generally 
consistent but that the reliability was variable across 
the frequency band (42). A recent rs-fMRI study using 
Pearson's CC analyses among raw time courses found that 
the correlation between the DMN and the 'anticorrelated' 
network of the resting-state network can vary over time 
(43). Thus, the functional connectivity of the resting- 
state brain network may be dynamically changed, yet may 
be consistent with the exception of a few RSIC pairs when 
functional connectivity is averaged within a certain period. 
The RSIC pairs showing the differences across time were 
more frequent in the CC analysis using Pearson's coefficient 
than in the PC analysis. Moreover, functional connectivity 
of the anterior DMN was variable across time in CC, but 
not PC. Furthermore, in CC, but not in PC, the upper part 
of the precuneus and inferior occipital gyri (occipital pole) 
showed differences in functional connectivity with multiple 
RSICs across time. This suggests that a possible indirect 
effect of functional connectivity among other RSICs on 
the functional connectivity of a certain RSIC pair can be 
controlled for using a PC analysis. A functional connectivity 
study using simulated fMRI data reported that a PC analysis 
was more reliable than a CC analysis (35). Consequently, the 
PC results are preferable in the present study. Interestingly, 
in both the PC and CC analysis, functional connectivity 
decreased across time between IC 12 (right angular and 
supramarginal gyri/intraparietal lobule) and IC 19 (bilateral 
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middle temporal gyri). IC 19 comprises bilateral V5 areas 
with connectivity to VI and V2 and the inferior and medial 
temporal gyrus. IC 19 is a region of the extrastriate visual 
cortex and is thought to play a role in the perception of 
motion and the guidance of some eye movements (44). In 
an fMRI study on a blind cohort, dorsal occipito-temporal 
regions were activated during the detection of auditory 
motion (45). The activated area in blindness was close to 
the maximum intensity point of IC 19 in the present study. 
Also, IC 12 was connected to the angular and supramarginal 
gyri and the visual cortex, consisted primarily of the 
parietal component of the right fronto-parietal network. 
The functional connectivity between IC 18 (posterior medial 
visual component) and IC 22 (right precentral/postcentral 
gyrus) was also variable in the PC analysis. Thus, there 
are some possible explanations for the variable functional 
connectivity of these visual information-process-related 
RSICs. Even if subjects have been asked to close their eyes, 
they could sometimes open their eyes, roll their eyeballs, or 
engage in visuospatial processing of scanner-induced noise 
during fMRI scanning. Their possible brain activity could 
modulate the functional connectivity of visual information- 
process-related RSICs. When subjects underwent fMRI 
scanning again, they might be less exploratory and more 
comfortable with their environment. If so, their functional 
connectivity under the relatively unfamiliar environment at 
day one could be reduced at 4 weeks. This may be a reason 
that the negative (IC 12-IC 19: bluish) or positive (IC 18- 
IC 22: yellowish) functional connectivity among RSICs 'at 
day 1' was not significant (whitish) at the 'after 4 weeks' 
time point in the present study (Fig. 2C). Our study showed 
no difference in functional connectivity across time after 
performing the correction for multiple comparisons. This 
should be not interpreted as implying that all functional 
connections are consistent. Rather, we argue that the 
correction for multiple comparisons should be applied in 
functional connectivity studies to reduce the possible false- 
positive error rate. 

In our study, none of the RSICs showed inconsistency 
in either functional connectivity or spatial characteristics. 
However, in a disordered brain, the intensity of the change 
could be different between functional connectivity and 
the activation intensity of brain regions in a network. 
Thus, to understand the brain as a dynamic network, we 
would recommend the investigation of the modulation of 
functional connectivity in brain networks as well as the 
change in activation of the brain regions. 
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We should note some limitations of the present study. 
The spatial map of the identified RSICs could be slightly 
different at every ICA analysis, and thus their functional 
connectivity could be variable as well. Our results could be 
confirmed by ICA at various dimensions in a further study, 
as the number of RSICs produced by ICA differs according 
to the magnitude of the dimension used for data reduction. 
Inter-subject variance in areas such as mental status during 
scan, education, or intelligence might cause different 
functional connectivity among RSICs. Thus, our results 
require confirmation with a dataset from a large sample, 
even though one study verified ICA repeatability using 42 
RSICs (46). The neuropsychological or clinical meaning 
of the functional connectivity of each IC pair should be 
explored in further studies. 

We found that most of the RSICs were reproducible across 
time, whereas some RSICs were variable in both their spatial 
characteristics and functional connectivity. Functional 
connectivity might be affected by the correlation analysis 
method applied. Our results suggested that researchers 
should consider both the variability of functional 
connectivity among the RSICs across time and the influence 
of the correlation analysis method. 
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